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Abstract 

New approach to classification of integrable hydrodynamic chains is established. 
Generating functions of conservation laws are classified by the method of hydro- 
dynamic reductions. N parametric family of explicit hydrodynamic reductions al- 
lows to reconstruct corresponding hydrodynamic chains. Plenty new hydrodynamic 
chains are found. 
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1 Introduction 

The first integrable hydro dynamic chain 

A'l = Al+^ + kA^~^Al, A; = 0,1,2,... (1) 

was derived in [1] . The integrabihty of this hydrodynamic chain was developed in a set of 
pubhcations ([13], [15], [20], [33]). The "integrability" means the existence of infinitely 
many conservation laws 

a,Hfc(AO, A\..., A^) = A\ A'+'), = 0, 1, 2... (2) 

and infinitely many commuting flows (see [20]) 

A^^ = [kA'^+^^-'d, + nd,A'^+''-^]-^, k,n,m = 0,l,2,..., 



where x = t^,t = t\ Hq = A^, Hi = A\ H2 = + {A^f, H3 = + 3A°A\ ... 
The first commuting flow is 

= ^ + (A; + 1)AM° + fcA'^-iAi, A; = 0,1,2,..., (3) 

where y = t^ . Eliminating the moments A^ and A^ from the equations 

A° = Al A\ = dM' + (^°)V2], Al = + {Ay] 

one can obtain the Khohlov-Zabolotzkaya equation 

Al = dMl-A'A% (4) 

which is the dispersionless limit of the KP equation. The compatibility condition dt{py) = 
dy{pt) of two generating functions of conservation laws for the Benney hydrodynamic chain 

(1) 

'p2 



Pt = d4^ + A']. (5) 



and for its first commuting flow (3) 



(6) 



yields (4). Thus, the integrability of this equation is equivalent the integrability of the 
Benney hydrodynamic chain (see [3], [14], [17], [18]). Substituting the series 



Ho Hi H2 , , 

P = X---^-^-... (7) 



in both above equations one can obtain (2) and the infinite series of conservation laws for 
(3) 

dyiik{A', A\ A'^) = d,Qk{A'>, A\ A''+'), k = 0,l, 2... 
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Later plenty integrable hydro dynamic chains 



n 



/c = 0,l,2,..., 




0, m> k + 1 
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■X 1 



n=0 



and corresponding 2+1 quasihnear equations (see, for instance, (4)) was found in [10], 
[11], [9], [12], [19], [23], [32]). Recently some of these integrable hydrodynamic chains 
were rediscovered (see [2]) and studied in [2], [4], [5], [22], [25], [24]). 

At this moment we have several tools allowing to find a complete classification of the 
integrable hydrodynamic chains (8) and more complicated hydrodynamic chains 



where and are some integers. 

1. The existence of an extra commuting flow (an external method), successfully used 
for the Egorov hydrodynamic chains (see [23]); 

2. The existence of an extra conservation law (an external method), successfully 
(an incomplete classification) used for the Kupershmidt Poisson bracket and for the (an 
incomplete classification) Kupershmidt-Manin Poisson bracket (see [19]); 

3. Vanishing of the Haantjes tensor (an internal method), successfully (unfinished 
classification) used for the integrable hydrodynamic chains (8) and successfully (complete 
classification) for the Kupershmidt-Manin Poisson bracket (see [9]); 

In this paper we establish the combined method based on three key tools, previously 
successfully utilized in the theory of 2+1 hydrodynamic type systems and 2+1 quasilinear 
equations (see [10], [11], [12], [32]). These are method of hydrodynamic reductions (see 
[15], [10]) and the method of pseudopotentials (see [10], [23], [32]) combined with the 
"concept" of the Gibbons equation (see [27]). 

The paper is organized in the following order. In the second section the concept of the 
Gibbons equation is introduced. In comparison with ([15]) hydrodynamic reductions are 
described in a conservative form. In the third section a new approach of the classification 
of integrable hydrodynamic chains via their generating functions of conservation laws is 
established. In the fourth section explicit hydrodynamic reductions associated with lo- 
cal Hamiltonian structures are discussed. In the fifth section the simple (but nontrivial) 
generalization of the Benney hydrodynamic chain is investigated. N parametric family of 
the Hamiltonian hydrodynamic reductions is found. The corresponding Riemann surface 
is constructed. Thus, series of conservation laws and commuting fiows can be found 
(see [27]), then infinitely many particular solutions can be obtained by the generalized 
hodograph method (see [31]). In the sixth section another integrable hydrodynamic chain 
is described, whose hydrodynamic reductions coincide with hydrodynamic reductions of 
the hydrodynamic chain described in the fifth section. One particular case of this hydro- 
dynamic chain is connected with 2+1 quasilinear equation determined by the "integrable" 
Lagrangian considered in [12]. In the seventh section three different approaches allowing 
to construct commuting flows are presented. The Miura type transformation is used for 
links between some well-known hydrodynamic chains and corresponding 2+1 quasilinear 






0, m > Mfe, 



(9) 



n=0 
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systems via generating functions of conservation laws. In the eight section we briefly 
mention the general case of generating functions of conservation laws connected with ar- 
bitrary integrable hydrodynamic chains. This general case is closely connected with the 
Hamiltonian hydrodynamic chains considered in [28]. 



2 The Gibbons equation 

The Benney hydrodynamic chain (1) is connected with the formal series 

AO 

\ = p+ — + — + — + ... 10 

p p-^ 



by the Gibbons equation (see details below) 



At - PK r, 

op 



The method of hydrodynamic reductions suggested in [15] (and developed in [10]) means 
the existence of infinitely many sub-systems (which are called the "hydrodynamic reductions'' 

rl=p\v)rl ^ = l,2,...,iV, (11) 

where all moments A^ are functions of Riemann invariants . This hydrodynamic type 
system must be consistent with the generating function of conservation laws (5). Then 
one can obtain 

d,p = (12) 
The compatibility conditions di{dkp) = dk{dip) yield the system in involution 

OiP = — r, OikA^ = 2— 7--, I k, (13) 

pi _ pk ^pi _ pky 

which we call the Gibbons-Tsarev system. 

Any hydrodynamic reduction (11) can be written via first moments A^ {k = 0, 1, A^— 
1), where all other moments A"^ [n = N, N + 1, ...) are functions of the first A^ moments 
A^ (see details in [15]). 

Also, any hydrodynamic reduction (11) can be written in the conservative form (see 
(5)) 



al = 5x(^ + ^°(a)), z = l,2,...,Ar, (14) 



where the function A^ satisfies the Gibbons-Tsarev system (written via field variables a'^) 



(15) 



/ i k\ 9ikA^ , / k j\ '^jkA^ , / j i\ ^ij^^ n • ^ ■ ^ / 



which is a consequence of the compatibihty conditions di{dkp) = dk{dip) 



diA'^ p- a"^ dkA^ p - a"^ diA'^dkA^ p- a 

where di = d/dd' and (cf. (12)) 

p — a'' y-^ — ^ p — a'^ J 

Example: The simplest hydrodynamic reduction is given by A^ = ^Ekd^. This is 
so-called "waterbag" reduction (see [16], [17]). Then (17) 

d,p = ^^(y^^-i 

p — a' V"^^ p — 
can be integrated. The equation of the Riemann surface 

x = P-Y,^kHp-a''), (18) 

where A is an integration factor, can be expanded (at the infinity A oo, p oo) in the 
formal series (10), where A'' = T.ei{a'f^'^ / {k + 1), if Se^ = 0. If Se^ ^ 0, then at first 
the parameter A in the equation of the Riemann surface must be re-scaled 

A-^efclnA=p - ^efcln(p - a!"). 

Then all other moments A^ are some non-homogeneous polynomials of Eej(a')'^. 

This is not an unique choice. For instance, any hydrodynamic reductions can be 
written in a conservative form in infinitely many ways. Let us mention just two simplest 
choices here (all other are discussed in [27]): 

1. dt = 4 (^^ + A°(a,b)^ , U, = d^{a'U) z = 1, 2, iV, (19) 

2. aj = d,(^^ + A'^, A° = 9.Ai(a) ^ = 1, 2, TV, 

For the first choice, the simplest hydrodynamic reduction is given by A^ = S6^; for the 
second choice A^ = T,b^. 

The Riemann invariants are most suitable coordinates to prove such general properties 
like the existence of infinitely many hydrodynamic reductions parameterized by arbi- 
trary functions of a single variable (see [15]), but they are very inconvenient for a search 
of particular and explicit hydrodynamic reductions. Vice versa, it is very difficult to de- 
rive the Gibbons-Tsarev system in the coordinates a'' (see (17)), but explicit particular 
hydrodynamic reductions can be found naturally in many cases (see a lot of examples in 
[25]). 
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The phenomenological algebro-geometric approach for integrabihty of symmetric hy- 
drodynamic type systems 

al = d^ip{a^,a'^,...,a^;p)\p=a' (20) 

was formulated in [27]. 

Statement 1: // the symmetric hydrodynamic type system (20) is integrahle, then 
this system has the generating function of conservation laws 

Pt = d^ijj{a^,a'^,...,a^;p). (21) 

Statement 2: Then some function A(a^, a^, a^; p) satisfies the Gibbons equation 

Xt-^K = ^[Pt-dM^;p)]- (22) 

We call the function A(a^, a^, a^; p) the equation of the Riemann surface. This 
function is a solution of the set of linear PDE's 

* du'' du^ dp 

where the matrix Af is given by 

The Gibbons equation (22) describes a deformation of the Riemann surface A(a;p). 
This equation has three distinguish features: 

1. if one fixes A = const (free parameter), then one obtains (21), 

2. if one fixes p = const (free parameter), then one obtains the kinetic equation (a 
collisionless Vlasov equation) written in so-called Lax form 

X _ M H} - - 

' dx dp dp dx ' 

where H = ip{a]p). 

3. if one choose coordinates, which are the Riemann invariants determined by the 
condition dX/dp = 0, then the corresponding hydrodynamic type system (20) can be 
written in the diagonal form (cf. (11)) 

rl,=i2\r)rl, z = 1, 2, iV, (23) 

where the characteristic velocities 



dp 



p=p'- 



can be found from the algebraic system dety4^'(u;p) = 0; then the corresponding values 
can be expressed via the Riemann invariants r'^. In this algebro-geometric construction 
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the Riemann invariants are the branch points r* = X\dx/dp=o of the Riemann surface 
(exactly as it is in the Whitham theory, see [7] and [18]). 

Example: The dispersionless limit of the vector NLS (see (19) and [33]) 

ai = d,(^^ + J2b^y bi = d.ia'b'), z = l,2,...,iV (24) 

is the first known hydrodynamic reduction of the Benney hydrodynamic chain (1) deter- 
mined by the moment decomposition A'^ = ^(a*)'^6*. In such case the formal series (10) 
reduces to the equation of the Riemann surface 

A = P + E-^- (25) 
The Zakharov reduction (24) written in the Riemann invariants (23) (see (11)) 

r\=p\v)rl ^ = l,2,...,2iV, 
where the characteristic velocities p*(r) can be found from [dX/dp = 0) 



E 



{p — a^y 



has infinitely many conservation laws, which can be obtained from (25) with the aid of 
the Biirmann-Lagrange series (see, for instance, [21] and some details in [27]). 

The compatibility conditions di{dkp) = dk{dip) must be valid identically for any sym- 
metric hydrodynamic type system (20), but a whole set (parameterized by arbitrary 
functions of a single variable) of such symmetric hydrodynamic type systems (for each 
fixed function ?/^(a;p)) is described by these compatibility conditions di{dkp) = dk{dip). 

Thus, the problem of the description of semi-Hamiltonian symmetric hydrodynamic 
type system is the problem of the classification of integrable hydrodynamic chains. 



3 The Gibbons— Tsarev system 

The main claim of this paper is that the classification of integrable hydrodynamic chains 
(9) is equivalent to the classification of the generating functions of conservation laws (21) 
by the method of hydrodynamic reductions (see [15], [10]). 

In this section we derive integrability conditions, which are a nonlinear PDE's sys- 
tem in involution. This extended Gibbons-Tsarev system generalizes the Gibbons-Tsarev 
system obtained in [15] and describes a set of generating functions of conservation laws 
together with their hydrodynamic reductions parameterized by arbitrary functions of 
a single variable. 

The simplest case is 

Pt = d^tp{u;p). (26) 
Let us first introduce the new notations 

/, = f{u,p\p) = , ^ f{u,p\p') = '^"^=;'' , t + k, (27) 

Yp\p=p^ Yp Tp\p=p^ Yp\p=p'' 
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,n (o, rr^i _ fikdpkfk - fkidp^fi + fjdpfk - fkdpfi - duifi - fk) , , . . 

Ji Jk 

Differentiating the generating function of conservation laws (26) witli respect to tlie Rie- 
mann invariants r* (see (23)), one obtains 

dip = fidiU. (29) 

U p = {k ^ i), tlien (29) reduces to 

V = UdiU. (30) 

The compatibihty condition dk{dip) = di{dkp) yields 

dikU = iPi^diudkU. (31) 

Thus, the nonlinear PDE's system (30), (31) is in involution iff the functions ip^^. do not 
depend on p and the compatibility conditions dj{dip^) = di{djp^) and dj{diku) = di{djku) 
are fulfilled identically. We call the system (30), (31) the Gibbons-Tsarev system (cf. 
(13)). 

However, in general case the functions (Pii^ depend on p. The compatibility conditions 
dkidip^) = d,{dkP') yields (cf. (31)) 

dikU = iPii^diudkU, (32) 

where (cf. (28)) 

¥>ik{'^,p) = <^ikiu,P\p'',P)\p=pi, iT^JT^k. 

Thus, the Gibbons-Tsarev system is determined by (30) and (32). 

Finally, one must check the compatibility conditions dj{diku) = di{djku) 

dui^ik-'Pjk)+'Piji'Pik-'Pjk)+fjidp^'Pik+fjkdpk(Pik+fjdpipik = fijdpj^jk+fikdpkcpjk+fidp^jk 

following from (31), and the compatibility conditions dj{diku) = di{djku) following from 
(32). Such over-determined nonlinear PDE system is said to be the extended Gibbons- 
Tsarev system. 

Let me emphasize that the extended Gibbons-Tsarev system is a system on the sole 
function ipiu^p) only. The general solution of this system yields a classification of inte- 
grable hydrodynamic chains. 

Remark: Any N hydrodynamic reduction (23) can be written in the conservative 
form (see (26)) 

aj = a,V'(w(a);aO, ^ = 1, 2, iV. (33) 

If the function ilj{u{a.);p) determines the integrable hydrodynamic type system (33), then 
the compatibility conditions di{dkp) = dk{dip) satisfy identically, where dk = d/da^ and 
(cf. (17)) 



diP 



(34) 



If the function ip{u;p) determines the integrable hydro dynamic chain, then the compat- 
ibihty conditions di{dkp) = dk{dip) describe its N component hydrodynamic reductions 
parameterized by arbitrary functions of a single variable. If the function ip{u]p) is 
unknown, then the compatibility conditions di{dkp) = dk{dip) and the compatibility con- 
ditions dj{diku) = di{djku) yield the extended Gibbons-Tsarev system. 

Let us finally emphasize: the Gibbons-Tsarev system describes component hydro- 
dynamic reductions parameterized by arbitrary functions of a single variable for any a 
priori given function ilj{u;p); the extended Gibbons-Tsarev system describes all possible 
functions ip{u]p) and their component hydrodynamic reductions parameterized by A^ 
arbitrary functions of a single variable. 



4 Explicit Hamiltonian hydrodynamic reductions 

Suppose the hydrodynamic reductions (14) of the Benney hydrodynamic chain (1) (si- 
multaneously, they are hydrodynamic reductions of the Khohlov-Zabolotzkaya equation 
(4)) are Hamiltonian (see details in [8]) 

where Si are arbitrary constants. Then the Hamiltonian density h = TjSkia^Y /Q + f{^)i 
where A = Ssfca*^ and /'(A) = A°(A). Substitution ^^(A) in (16) yields the choice 
A'' (A) = A only (up to an insufficient constant factor). This is exactly so-called "wa- 
terbag" hydrodynamic reduction (see (18)). 

Thus, the main claim of this section is that the Hamiltonian hydrodynamic reductions 
(cf. (33)) 

a\ = ^^^lJ{u■, a') 

for any hydrodynamic chain determined by the the generating function of conservation 
laws (see (22) and (26)) 

pt = d^i/j{u;p), 

one can seek in the form 

al = dJS,^ + jyik^], (35) 



where Si and 7^ are some constants (see examples in [27]). 

The second step is a reconstruction of the Riemann surface determined by the equation 
A(u;p) by virtue of (35) and a computation of the moments = E/^(a"'), where the 
functions /^(a") from the expansion of the equation of the Riemann surface at infinity 
{\ ^ oo,p ^ 00). Couple of examples are considered below. 



5 Generalized Benney hydrodynamic chain 

In this section we restrict our consideration on the first simple and very important par- 
ticular case (cf. (5)) 

Pt = d,[U{p)+u]. (36) 
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Theorem: component hydrodynamic reductions (23) 

r\ = u'ipyi 

are described by solutions of the Gibbons-Tsarev system (cf. (13)) 
. _ 2aU'{p')U'{p>^) + f3[U'{f) + U'jp'^)] +27 ^ 



(37) 



U'{p') - U'{p) ' U'{p^) - U'{p^ 



where the function U{p) is a solution of the second order ODE with respect to the inde- 
pendent variable p 

U"{p) = aU'\p) + pU\p) + 7. (38) 

The compatibihty conditions dj{dikA^) = di{djkA^) are identicaUy satisfied for any 

constants a, (5, 7. In general case (a 7^ 0) the function U{p) can be written in the 
parametric form only 

U = ——^ -[giln(g-gi) - ln(g - ^2)], P= —r-^ -[In(g-gi) -ln(g-g2)]- 

a{qi - q2) a{qi - ^2) 

If a = 0, but /9 7^ 0, then without lost of generality U = e^; if /3 = 0, then U = if 
a 7^ 0, but P = and 7 = 0, then U = Inp; if a 7^ 0, but /? = 0, then U = Insinhp (or 
f/ = In cosh p) . 

The first step in the reconstruction of integrable hydrodynamic chains is 
an extracting of any explicit hydrodynamic reduction (see (33)) 

al = d,[U{a') + u{s.)]. (39) 

Then the compatibility conditions di{dkp) = dk{dip) yields (cf. (15), (16)) (38), two 
distinct index relationship 

[U\a') - U'ia')]{d^kU + ad.udku) = duudi ( J] - d^udu (j] , (40) 

and three distinct index relationship 

[f/'(a') - U'{a^)]djudikU + [U\a^) - U'{a^)]diudjkU + [U'{a^) - U' {a')]dkud,jU = 0, 

where (see (34)) 



d.p= ^ ^ + 1] (41) 

Following the previous section we are looking for the Hamiltonian hydrodynamic re- 
ductions (39) 

4 = d,[U{a')+u{^)] = jd,^^. (42) 
10 
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Theorem: The Hamiltonian hydrodynamic reductions (42) are determined by the 
Hamiltonian density h = Se^ jU{a^)dal' + f{A), where A = Se^a'^, /'(A) = u{A) and 
u{A) = A, z/ a = 0; n(A) = ln(A), zf a ^ 0. 

Proof: can be obtained by the substitution (42) in (40). 

The equation of the Riemann surface (see (41)) can be found in quadratures 

n'(A) ^ ^ [/'(a'^) - t/'(p) ^ ' 

for this Hamihonian hydrodynamic reduction. 

Let us introduce the moments = Se^ J U' {a"')da'^, then the hydrodynamic type 
system (42) can be rewritten as the Hamihonian hydrodynamic chain 

fc+2 

Al = J2F!:iA)A:, k = 0,1,2,... 

n=0 

determined by the Hamihonian density h(A°, A^) = A^ + /(A°) (where /(A°) = {A^y/2, 
if a = and f{A^) = A^lnA'^) and by the Dorfman Poisson bracket (see [6]) 

{A\A'} = J2^nS'{x-x'), 

{A^ A"} = [A;(aA^'+"+^+/3A'=+"+7A'=+"-i)9^.+n9^(av4'=+"+^+/3A'=+"+7v4'=+"-i)](5(x-,x'). 
For instance, this hydrodynamic chain is 

A° = d4aA^ + PA' + (7 + 5^£n) A'] 

A\ = {aA''+^ + pA^+' + -iA^)^ + k{aA^^' + pA^ + 7^'=-^)^°, k = l,2, 
if the Hamihonian density is h(A°, A^) = + {A^Y/2. 

6 The Gibbons— Tsarev system and hydrodynamic chains 

In this section we consider the second simple and very important particular case 

Pt = d^[V{p)v]. (44) 

Then TV component hydrodynamic reductions (23) 

r\ = vV'{p')ri 

are compatible with the above generating function of conservation laws iff (see (29)) 

^'P = ^^p \w-npy 

It is easy to see that under re-scaling 

^ = dp, V'ip) = U'ip), 
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one get the main formula from the previous example (see (37)), where u = In v. The 
corresponding equation is (see (38)) 



VV" = aV' + (3V' + 7, 



(45) 



where 

V{j)) = exp U (p), dp = exp U {p)dp. (46) 
In general case (a ^ 0) the function V{p) can be written in the parametric form only 



V{p) = {q — Qi) "'■n-12) (^q — "(n-n) . 



1 f q-j 92 

p=— / (g — gi)°(''i-«2) (g — ^2) "("^i-'a) dq. 
a J 



The corresponding degenerate cases (see the above example) are V{p) = p'', where k is 
an arbitrary constant, V{p) = e^, V{p) = coshp (or V{p) = sinhp) and two others given 
in the implicit form 



p 



dV 



P 



dV 



Since the equations describing hydrodynamic reductions of these both examples (36) 
and (44) are equivalent, then one can recalculate any hydrodynamic reduction of the 
one example to a corresponding hydrodynamic reduction of the second example (the 
Hamiltonian reductions of the generating function of conservation laws (44) are found 
in [27]). The relationship (46) is nontrivial. For instance, U{p) = p^/2 for the Benney 
hydrodynamic chain (see (5)), but V{p) cannot be expressed via known elementary or 
special functions. Nevertheless, hydrodynamic reductions of the Benney hydrodynamic 
chain can be recalculated. One can compare the formulas (34) for (36) and (44) 



dp V{p)dv/dd 

da^ ~ V'{a^) - V'{p) 

dp du / dd 

dd ^ U'{d) - U'{p) 



V 



y- 



V{a^)dv/da^ 



EC I 
TT'(rk 



V'{a^) -V'{p)_ 
du/dc'^ 



U'{&) - U'{p) 



They are equivalent under the transformation (46) and under the same transformation 
for the field variables ^ 

V{a') = exp U{d), da' = exp U{d)dc\ 

The equation of the Riemann surface (43) also can be recalculated 



dX = (\/(p))'"-^exp[/5 



dp f v{A)dp 



V{py \ v'{A) 



Ek V{a^)dp - V{p)da}' 



where A = Se^ / da^/V{a^), v{A) = exp A, if a = 0; v{A) = A, if a 7^ 0. 

Example: The Lagrangian quasilinear equation associated with an elliptic curve. 
The Lagrangian 

ZxZyZtdxdydt 
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creates the Euler-Lagrange equation 



yt 



(47) 



which is an integrable 2+1 quasihnear equation (see [12]). The pseudopotentials (cf. (4) 
(5) and (6); see also [10], [32]) are 



2p(a) 



(48) 



where C,{a) and p'{(t) are Weiershtrass eUiptic functions (C'(cr) = ~p(c"), p'^(o") = 4p^(cr) + 
e^). Introducing the new functions a = Zx-, h = Zy, c = Zt, (47) can be written in the form 



bt Cyj 



ctty + bat + ctbt = 0. 



Correspondingly, (48) can be written as the couple of generating functions of conservation 
laws 



Py = dx 



a ^ 2p{a) 



a 



where p = Sx and a can be found from ({a) = p/a. 
The reciprocal transformation (see [29], [30]) 

dz = adx + bdy + cdt, dt = dt, 

reduces the above Lagrangian to 



2p{a) 



dy = dy 



(49) 



dzdydi 



and (49) to the couple of generating functions of conservation laws 



Py = 



2p{a) 



Pi = dz 



<pip)-^ 
2p[a) 



(50) 



where p = C(^)- However, these generating functions of conservation laws are considered 
in this section (see (44)). Indeed, the equation (45) 



V 



in 2 



3 



\ dp 



dV 
+ 9— + 6 



has the general solution given in the parametric form 

p'{a)±e 



V{p) 



2p{a) 



P = C(^) 
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7 Commuting flows and 2+1 quasilinear systems 

In this section three different approaches in construction of commuting ffows for any given 
hydro dynamic chain are presented. 

1. Suppose all generating functions of conservation laws (26) are enumerated. Then 
generating functions of conservation laws for commuting flows one should seek in the 
forms 

Pti = d^ip^{u^;p), pt2 = d^ip2iu^,u'^;p), ... (51) 

The compatibility conditions dt{pt^) = dti{pt), dti^pt^) = dt^ipt), ... allow to reconstruct 
functions ipj^ in quadratures (the corresponding example is given by the Benney hydrody- 
namic chain (1) and the Khohlov-Zabolotzkaya equation (4); see (5) and (6)). 

Example: Let us consider the couple of commuting generating functions of conserva- 
tion laws (cf. (50)) 

Py = [V {p)b], pt = [W{p)c]. 
The compatibility condition dt{py) = dy{pt) yields (45) 

VV" = (l + ^)V'" + Yl-6-2^]V' + -{(35-a^), 
a a a 

WW" = {l + -)W'' + [a- (3-2^]W' + -{f36-a-i), 

7 7 7 

where 2+1 quasilinear system is (a,/3, 7,5 are arbitrary constants) 

ht = ahc^ + (3cb^, Cy = -fcb^ + 6bc^. 

2. Let us replace p{Q and t{Q in the first equation (51) 

drioPi\)=d,i;MC),pW)- (52) 

Definition: The equation (52) is called the generating function of conservation laws 
and commuting flows. 

All generating functions of conservation laws (51) can be obtained by the expansion 
of (52) in the series according the Biirmann-Lagrange series of p{C), where dr{(^) is the 
formal series, whose coefficients d^k enumerate different commuting flows. 

Example: The Benney hydrodynamic chain has the generating function of conser- 
vation laws (5) and the generating function of conservation laws and commuting flows 
is 

drioPW=d,\n[p{\)-p{C)], (53) 

where ^^(^) = dxPiO- Substituting the expansion (7) for p{\) in the above formula, one 
can obtain the generating function of commuting flows in the form 



14 



Substituting the expansion (7) for p{X) in the above formula, one can obtain the generating 
functions of conservation laws in the form (see (5) and (6)) 



3. Suppose for any given function ■ip{u] p) (see (26)) we already know the corresponding 
Hamiltonian hydrodynamic reductions 



da'' 

where g^'^ is a constant non-degenerate symmetric matrix. Then one can seek the higher 
conservation law density h2 for this hydrodynamic type system. This conservation law 
density determines the higher commuting flow 

The generating function of conservation laws for this hydrodynamic type system (as well 
as for the corresponding hydrodynamic chain) can be found by the replacement a* ^ p 
(see details in [27]). The generating functions of commuting flows is given by 

«r(C) - [9 ) ■ 

Let us replace dr[() dt-i and, correspondingly, p{() then the generating 

function of conservation laws (53) can be written in the form 

p. = dMe^ + A-'), (54) 

where the generating function of the Miura type transformations (see [26]) is given by 
p = ln{p — A^^). This generating function (54) determines the continuum limit of the 
discrete KP hierarchy (see [17]) 

Bl = 3^+^ + kB^Bl, A; = 0,1,2,..., 

where z = t~^ and = . 

Remark: All moments A'^ are connected with the moments B'^ by the Miura type 
transformations 



A^ = A^B'^^B^), A' = A\B'>,B\B^), A^ = A\B'',B\B^,B'^),..., 



which can be obtained by a substitution the generating function of the Miura type trans- 
formations p = A~^ -|- expp in (10) and a comparison with (see [17]) 

A = + 5° + B^e~^ + E^e'^^ + ... 
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Let us substitute the expansion d^-i^Q = dt-i + (dt~2 + ... and p{C) — > + CA~'^ + ... 
in (53). Then the compatibihty condition dt-2{dtip) = dti{dt-2p) yields the 2+1 shallow 
water system (see [32]) 

d,-2A' = d,A-\ d,^A~' = + ' = d,{A~'A~'), (55) 

where 

df-2p = dx—r-, • 

The compatibility condition dt-2{dt-ip) = dt-i{dt-2p) yields the famous Boyer-Finley 
equation 

dt-iA~^ =dAr^A-\ dt-iA"^ = dt-2A~\ (56) 

In the both case (55) and (56) TV component hydrodynamic reductions (11) and (see 
(23)) 

1 4-2 

p'-A-^ {p'-A-^y " 

are compatible (i.e. (9ti(rJ_2) = 9j-2(r*i); see [31], [10]) 

dk ( „i \-i ] dk 



iff the functions and A^ satisfy to the Gibbons-Tsarev system (13), where 



It is not easy to verify, if did not take into account, that all these 2+1 quasilinear systems 
(4), (55) and (56) are members of the unique Benney hydrodynamic lattice (see [25]). 



8 General case 

The main statement of this paper is that all integrable hydrodynamic chains described 
by the generating functions of conservation laws (21) can be split into sub-classes accord- 
ing to the number M of the functions u™. The first such a case (26) was considered in 
the previous section (M = 1). The second case (M = 2) is 

Pt = dxij{u,v;p). 

Then component hydrodynamic reductions (23) are compatible with this generating 
function if the compatibility conditions di{dkp) = dk{dip) are fulfilled, where 
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To complete this computation, we need some relationship between u and v (except trivial 
link v{u)). Without lost of generality we can choose these field variables u and f as a 
conservation law density and a flux, respectively: 



Ut = v^. 

Since the hydrodynamic type system (23) has this conservation law, then 

diV = ^p\p=pidiU, 

and (cf. (27), (29)) 

However, in such a case the Gibbons-Tsarev system can be derives from the extended 
compatibility conditions di{dkp) = dk{dip) and di{dkv) = dk{div). 



9 Conclusion and outlook 

In this paper a new look on a classification of the integrable hydrodynamic chains is 
presented. The main object which should be under an investigation is the generating 
function of conservation laws 

Pt = d^ij{u;p), 

where all distinct cases are separated by the number M of independent functions u™. The 
simplest case (26) is considered in details. Each such case has infinitely many sub-cases 
enumerated by the number K of independent functions in the commuting generating 
functions of conservation laws 

Py = d^ip{v;p). 

Thus, all integrable hydrodynamic chains can be split into sub-classes by virtue of the 
two numbers M and K only. 

Suppose all these functions iplu^p) are found. Then the Gibbons-Tsarev system de- 
scribing N component hydrodynamic reductions for every function ?/^(u;p) can be derived 
automatically in the Riemann invariants and in the field variables (which are conserva- 
tion law densities) simultaneously. Then corresponding hydrodynamic type systems with 
a local Hamiltonian structure (as well as with a priori prescribed any nonlocal Hamilto- 
nian structure) can be extracted. The equation of the Riemann surface X{u;p) can be 
found in quadratures for the hydrodynamic reductions, whose characteristic velocities are 
invariant with respect to any Lie group symmetry. Asymptotic of the equation of the Rie- 
mann surface at the vicinity of any singular point (usually A —>■ oo,p —>■ oo) yields explicit 
expressions of corresponding integrable hydrodynamic chains. The "fiat" hydrodynamic 
reductions can be used for a construction of a large class of particular solutions for these 
hydrodynamic chains and related 2+1 quasilinear equations. 
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